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Abstract 

We introduce a shell ("GOY") model for turbulent binary fluids. The variation in the 
concentration between the two fluids acts as an active scalar leading to a redefined conservation 
law for the energy, which is incorporated into the model together with a conservation law for 
the scalar. The model is studied numerically at very high values of the Prandtl and Reynolds 
numbers and we investigate the properties close to the critical point of the miscibility gap 
where the diffusivity vanishes. A peak develops in the spectrum of the scalar, showing that 
a strongly turbulent flow leads to an increase in the mixing time. The peak is, however, not 
very pronounced. The mixing time diverges with the Prandtl number as a power law with an 
exponent « 0.9. The continuum limit of the shell equations leads to a set of equations which can 
be solved by a scaling ansatz, consistent with an exact scaling of the Navier-Stokes equations in 
the inertial range. In this case a weak peak also persists for a certain time in the spectrum of 
the scalar. Exact analytic solutions of the continuous shell equations are derived in the inertial 
range. Starting with fluids at rest, from an initial variation of the concentration difference, 
one can provoke a "spontaneous" generation of a velocity field, analogous to MHD in the early 
universe. 
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1 Introduction: Turbulent fluid mixtures. 



Binary fluid mixtures provide a beautiful example of physical systems where it is possible to study 
the behavior in the limit of exceedingly large values of the Prandtl number. This is the case for 
miscible binary mixtures just above the consolute temperature T c , at which the diffusivity vanishes 
as a power law Q 

D(T)^D (^-^Y (1) 



with an exponent in the range 

0.63 -0.74 . (2) 

Since D(T) — > at the critical point, the corresponding Prandtl number Pr = u/D — > oo and values 
as high as 10 6 can be obtained experimentally M. Goldburg and coworkers studied turbulent binary 
mixtures experimentally using light scattering techniques || ^, f§, ||. One can measure both the 
variation of the mixing times and the growth of domains using these techniques and we return to 
a discussion of the experimental results later. 

Our motivation for the present work is that we are able to formulate the theoretical equations 
behind binary fluid mixtures in terms of shell models where the known invariant quantities are 
conserved. In this way one can investigate the mixture at much higher values of the Reynolds and 
Prandtl numbers, than is possible by standard numerical simulations. For instance, one can then 
study the scaling behavior of the mixing time at very high values of the Prandtl number, a regime 
which might be experimentally accessible. Also, the shell model is known to include intermittency 
effects, which have not been treated in previous theoretical works on binary fluid mixtures. 

With the two fluids labelled A and B, respectively, a scalar field is defined as ||, |7j 

i>{r,t) = ( PA (r,t)-p B (r,t))/ Po (3) 

where pa(*", t) and p_B(r, t) are the mass densities of the two fluids and po is the mean mass density. 
In a phase plane determined by the temperature T versus the average of the scalar < ip >, there 
exists a "miscibility-gap" separating the miscible phase from the immiscible. Along this separating 
curve, the effective diffusivity D vanishes, because it separates a regime where the effective diffu- 
sivity is positive (i.e the miscible case) from a regime where it is negative (the immiscible case). 
Close to the 50-50 % concentrations of the two fluids, one finds in equilibrium 

<VM)> =o . (4) 

For this case, as the critical point is approached from the miscible phase, the scalar is supposed to 
be "active" and influence the velocity equation of the Navier-Stokes equations quite substantially. 
The corresponding equations of motion were derived almost two decades ago by Siggia, Halperin 
and Hohenberg (6|, ^] and later on elaborated quite a lot by Ruiz and Nelson [||, ||] 

^ + (u • V)^ = .DVV (5) 
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+ (u • V)u = Vp' - aWVV + ^V 2 u + f (6) 

at po 

V • u = (7) 

Here v is the kinematic viscosity; f is the forcing and several terms involving ip have been incorpo- 
rated into an effective pressure p' [||, [7| . The term with coefficient a represents the "active" part of 
the scalar. This term acts like a force of the form fiAB^tp where \xab = — aV 2 ip plays the role of 
a local chemical potential difference between the A and B component of the mixture H-|J. The 
coefficient a has the dimensions of the square of a transport coefficient and has been estimated to 
be of the order a ~ v 2 in reference 

The equations of motion and (^) allow two quadratic invariants in the absence of diffusivity, 
viscosity and forcing, i.e. in the limit D = v = 0, f = 0. The first is the squared integral of the 
concentration fluctuations 

C tot = \ fdr^(r,t)) 2 (8) 
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and the second is the total energy with a term relating to the active influence of the scalar 

Etot = \J dr{\ u(r, t) \ 2 +a | V%b(r, t) | 2 ) . (9) 

In the case of a passive scalar, i.e. when a = 0, one expects the energy spectrum E(k) and the 
spectrum of the scalar C{k) to have the usual behavior 

E(k) ~ AT 5/3 ~ 5 , C(k) ~ k- 5/3 ^ (10) 



where -5/3 is the Kolmogorov exponent (and Obukhov-Corrsin [11| exponent for the scalar) 
and 5 and 7 are intermittency corrections in the two cases, respectively [12]. Ruiz and Nelson [R 0] 
also discuss the possibility of internal wave-like excitations, in the case of large values of a, similar 
to linear wave excitations in MHD, which may change the spectrum to different scaling behavior 



as predicted by Iroshnikov and Kraichnan 13, 14]. We do not discuss this phenomenon here but 
reserve it for a forthcoming publication. 

An interesting feature of the equations of motion (||) and @ should be noticed: If the initial 
velocity field u vanishes for i=0, then if 

/ and V 2 ^ ± 0, and/or Vp' / for t = 0, (11) 

it follows from Eq. @ that a finite velocity field will appear. For small times it is given by 

u(x, t) = - (—Vp'{x, 0) + aVV>(x, 0)VV(x, 0) ] t + 0(t 2 ). (12) 
\Po J 

Here the pressure term should be such as to respect eq. (0). Eq. ( |i~2|) means that if initially the 
liquids are at rest, and experimental initial conditions respecting ([ll]) are established, then one 
should see that the liquids start to move "spontaneously". This is a very clean effect of a non- 
vanishing "transport coefficient" a, and it works irrespectively of the magnitude of the diffusion 
coefficient D. This will be further discussed in section 7. 
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The paper is organized as follows. In section 2 we derive the shell model for turbulent binary 
fluid mixtures and discuss the corresponding conservation law, the value of the coupling constants, 
etc. In Section 3 the numerical results obtained from integrating the model are presented. In 
particular we discuss the appearance of a peak in the spectrum of the scalar. Section 4 contains 
the theoretical predictions of Ruiz and Nelson for the mixing times and the corresponding results 
from the shell model. In section 5 we present the continuum version of the "GOY" model and exact 
analytic solutions in the inertial range, based on a scaling ansatz. In section 6 the corresponding 
continuum equations for the binary mixture model are derived and in section 7 the numerical results 
from integrating these equations are presented together with a comparison with the results from 
the discrete equations. Finally, section 8 offers concluding remarks. 



2 A shell model for binary mixtures. 

Since the binary mixtures are particularly interesting to investigate in the critical regime where 
Pr — > oo and as we are concerned with the case of a strongly turbulent mixture (large values of 
Re), it is our goal to formulate an approximate scheme for Eqs. (||) and (JgJ) in which this limit is 
accessible. Shell models in Fourier space fulfill these requirements. They have been introduced by 



Obukhov [15 1, Gledzer fllql , Desnyansky and Novikov [17]. The key idea is to mimic the Navier- 
Stokes equations by a dynamical system with N variables ui,U2, njy, each of which representing 
the typical magnitude of the velocity field on a certain length scale. The Fourier space is divided 
in N shells and each shell consists of the set of wavevectors k such that k^r 11 < | k | < kor n+1 . 
The variable u n is the velocity difference over a length ~ k~ l so that there is only one degree of 



freedom per shell [12|. Also models with a large number of degrees of freedom have been introduced 
and analysed |[^, The most studied model is the "GOY" model introduced by Ohkitani and 
Yamada [^0] which was found to be intermittent by Jensen, Paladin and Vulpiani [21| and studied 



extensively in many other contexts [23, 24, 25, £7], ^8], This model uses af complex set of 
variables and has the same type of quadratic non-linearities and the same symmetries as the 3d 
Navier-Stokes equations. We shall apply the same approach in this paper by expanding the "GOY" 
model to include the term of the active scalar and at the same time conserve the two different 
quadratic non-linearities (|8|) and (ph. 



Firstly, we write down the shell model for the scalar equation (|50 M, 22]. Using a complex field 
ip n associated to shell n, the equation becomes 

d 
•~dl 



(37 + Dkfy ^ n — i [e n k n (ll*_ 1 V'n+i — ^n+lVC-l) + 9n &n-l ( - U^_ 2 V'n-l + ^n-lV'n^) 



+ h n k n+1 « +1 ^+2 + < +2 VCk) ] (13) 

where D is the molecular diffusion. The physical time scale of the model is determined by the the 
constant k\ = rko, which represents the inverse scale of the largest eddy, and the related velocity 
\ui\. This time scale is therefore the corresponding eddy-turn-over time In the following, 

"time unit" thus means this time scale. 
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The coefficients of the advective terms follow from demanding the conservation of ^2 n \ip n \ 2 
when the diffusivity is vanishing D = 0. A possible choice is: 

1 1,1 

e n = - 9n = — K = - (14) 
r r r 

with 

ei = e N = gi = g 2 = h N -\ = h N = (15) 

The shell model equations for (g) will consist of two contributions; the first part is the usual 
"GOY" shell model p(J (with coefficients a n , b n , c n ) and the second part is the shell expression for 
the "active" term in the velocity equation @ 

(— + vk 2 ) u n = i k n [a n Un+l u n+2 + b n u* n _iU* n+l + C n M*_ 1 lt*_ 2 ] 

+ tnC-lC-2] + /<W (16) 

with n = 1, • • • N, k n = r n k$ and boundary conditions 

h = 6jv = ci = C2 = ajv-i = oat = si = sat = ti = t 2 = rjv_i = = (17) 
In order to ensure the conservation of the quadratic quantity 

Etot = E (l U « I 2 +ak n I I 2 ) (18) 



in the limit without viscosity, diffusivity and forcing, v = D = f = 0, one multiplies Eq.(16) by u n 
and multiplies Eq.(|l3|) by ^ n and then balances the terms. The non- linear terms in u n lead to the 
usual constraints of the "GOY" model [ fL^ ] 

6 - 1 Z 6 

For the second part of Eq.(|T6|) one balances the terms by the corresponding terms in the scalar 
equation. With the choice of the parameters (14) we then obtain the following conditions for the 
coefficients 

t n = r- 4 -r~ 6 . (20) 



One observes that when the coupling constants for the scalar equation are given (14), then the 
coefficients of the "active" terms are fixed. As these active terms are proportional to k^, we expect 
the effects of this term to show up at the end of the spectrum for large values of k n . 



3 Results from the shell model. 

This section contains some of the numerical results obtained from integrations of the shell model 
for the binary mixture, derived above. In the simulations we use the standard separation between 
the shells, r = 2, such that k n = 2 n /co- We apply the "symmetric choice" of the "GOY" parameters, 
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5 = 1/2. In this case it is known that for the "GOY" model alone, the second quadratic invariant 
assumes the symmetry of a helicity p3fl , and for that case the model is strongly intermittent and 
gives results in good agreement with experiments. We study the model with N = 14 and 19 shells, 
ko = 2~ 4 , and the strength of the forcing term in ( |i~6| ) is / = 0.005 • (1 + i). As argued by Ruiz 
and Nelson [j|, |9| the coupling constant of the active term a, has the dimensions of a square of a 
transport coefficient and is in the order of magnitude 

a ~ v 2 . (21) 

The spectrum of the scalar for the shell variables is defined as 

C(k n ) = <\ij n \ 2 >/k n . (22) 

The brackets stand for averages over initial conditions and time. Similarly, the energy spectrum is 
defined as 

E(k n ) = <\ u n | 2 > IK . (23) 

Since the scalar equations (|E|),(|l3"|) are not forced, then C(k n ) — ► in the long time limit. Nev- 
ertheless, it is possible to obtain intermediate averages over shorter times. This is in contrast to 
the case of the velocity spectrum E{k n ) where the mean exists for t — > oo as the velocity equations 
(PI) ; fll^D are forced. Firstly, we present results when the value of the viscosity is v = 10 -4 , meaning 
that Re ~ 10 4 . For this value of the viscosity, a shell model with iV=14 shells is employed. In 
order to observe the differences between a passive and an active scalar, we first consider the case 
a = 0. Fig. [l]a shows the corresponding spectrum C(k n ) on logarithmic scales in the case Pr = 1. 
The spectrum follows quite closely the Obukhov law pT| . 

Next the Prandtl number is increased to Pr = 10 3 . The spectrum, shown in Fig. [l|b, is changed 
and scales for high value of k n according to the Batchelor law C(k n ) ~ A;" 1 , which means that with 
the increased value of the Prandtl number a viscous-convective regime is observed, as expected ||. 
When the scalar becomes active, i.e. a ^ 0, by introducing the coupling term (p|), the spectrum 
changes completely and a peak develops at the upper end of the spectrum as indicated on the 
evolution series, Fig. 0c. The series is initiated in a state where Vn is concentrated at a low 
k n value. This corresponds to a large scale disturbance, for instance where the fluids A and B 
are completely separated. After some time, a peak (or, perhaps more appropriate, a "shoulder") 
develops at large values of k n as indicated in the figure (each curve is averaged over 200 time 
units). After further time the peak disappears and the spectrum ends up in C(k n ) ~ fc^ 1 . The 
final spectrum in Fig. 0c represents again an almost stationary situation and does not change 
significantly during long time. Let us note in passing that time-averaging give results similar to 
ensemble- averaging, since the dynamics in the phase space of the u n - and ^-field is strongly 



chaotic as determined by positive Lyapunov exponents [12]. 

This peak was predicted by Ruiz and Nelson and also seen in numerical simulations using 
Markorvian closure equations ||. In fact, in those simulations the peak appears much more pro- 
nounced than compared to our results. We believe the reason is that the strongly intermittent 
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Figure 1: Results from numerical integrations of the shell model and (|l6|). The parameters 
are: N=14, ko = 2~ 4 , v = 10 -4 . The spectrum in (a) is averaged over 4000 time units; the other 
over 200 time units. In the case (a) the spectrum C(k n ) has D =10 -4 , Pr =1, and a = 0. The 
dashed line in (a) has slope -5/3. In the case (b) D=10 -7 , Pr=10 3 , and a=0, and the dashed 
lines have slopes -1 and -5/3, respectively, (c): C(k n ) for D = 10~ 7 and a = 10~ 8 . The full 
curve shows the initial condition. Time progresses from the uppermost spectrum to the lower. The 
interval between the curves are of the order ~2000 time units. In (d) we show the spectrum at one 
particular time. The dashed line has slope -5/3. 
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motion could influence the dynamics in a way that the peak becomes less strong. Fig. |l|d shows 
one of the spectra in the series indicating the peak more clearly. Note, that for low & n -values, the 
spectrum is still close to the — | law as indicated by the dashed line. The peak indicates, that 
two miscible fluids close to the critical point mix very slowly at the small scales when the fluid is 
strongly turbulent. The time for which the peak persists is strongly dependent on the Prandtl and 
Reynolds numbers. This time is called the mixing time |]]. 

Fig. |2|a shows the development of the spectrum C(k n ) for a lower value of the viscosity, v = 1CP 6 




-10 12 3 4 -10 12 3 4 5 

L °gl0 k n L°gl0 k n 



Figure 2: To the left, spectra C{k n ) obtained from numerical integrations of the shell model, with 
parameters N = 19, k^= 2 -4 , v = 10 -6 , D = 10 -13 , a = 10~ 12 . The different spectra are averaged 
over 4 time units and the distance between the spectra are ~ 50 time units. Time progresses from 
the left to the right. The right hand figure shows one particular spectrum. The dashed line has 
slope -5/3. 

corresponding to Re ~ 10 6 . In this case D = 10" 13 leading to Pr = 10 7 . A gain, starting from an 
initial condition concentrated on the small /c n -values, one observes the occurrence of the peak at 
large /c n -values. Fig. shows one of the spectra and the dashed line corresponds to C(k n ) ~ k n 5 ^ 3 . 
In order to get a more clear picture of the peak, we plot <| ip n |> versus k n on logarithmic scales. 

— 1/3 

For a usual Obukhov spectrum one should find <| if) n |>~ k n , whereas for the viscous-convective 
regime <| ip n |>~ const, and the peak will therefore appear more pronounced. Fig. || shows the 
time development of the value of the scalar, where each curve is averaged over 800 time units 
(the time intervals between the curves are ~8000 time units). One observes that ip n is almost 
constant over the leftmost regime of the spectrum, corresponding to C(k n ) ~ A;" 1 , whereas the 
peak is situated at the rightmost part of the spectrum. The peak slowly decreases in intensity 
(at a rate determined by the mixing time) but for high values of the Prandtl number this decrease 
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Figure 3: A plot of <| ip n |> versus k n on logarithmic scales. The parameters are ./V = 19, ko = 
2 -4 , v = 10~ 6 , D = 1CP 12 , and a = 10 -12 . Each curve is averaged over 800 time units and the time 
intervals between the curves are ~8000 time units. Time progresses from the uppermost curve to 
the lower. 

occurs extremely slowly: the peak in Fig. [3| diminishes significantly only after about 10 6 time units. 
Simultaneously, the field as a whole slowly vanishes, ip n — ► 0. This is of course related to the fact 
that the passive scalar equation (|l^) is not forced and the input to the motion of tp n is only driven 
by the advective term which includes u n . Nevertheless, as Pr — > 0, ip n — ► slowly, which indeed 
is reflected in the corresponding value of the mixing time. The presence of a peak appears to be 
independent of the initial state; one can either, as in Figs.||j2|, apply an initial disturbance which 
is concentrated at the small k n values, or choose states of ip n and u n which are solutions to the 
passive scalar equations, i.e. in which a = 0. In all cases the results is a peak in ip n at large k n , 
so an initial perturbation at the large scales in not a necessity in order to observe the enhanced 
"delay" in the mixing of the two fluids. 

In comparison with previous work ||| and H, it should be emphasized that in our case the peak 
is orders of magnitude smaller in height. In Ref. ||, the peak (after a short time) has a height of 
order some decades. In our case, as seen from the various figures, the height is only of order half 
a decade at most, and in many cases it is much less, which could make experimental observations 
difficult, as indeed appears to be the case || 

The corresponding velocity spectrum does not show any sign of a peak, Fig. |I[ On the contrary, 
the cut-off at high wavenumbers in the spectrum in fact seems to move to shorter wavenumbers 
(larger scales) as a results of the active coupling. The oscillations in the spectrum is on the other 
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Figure 4: The velocity spectrum | u n | versus k n for the same parameters as in Fig. ||. The 
spectrum is averaged over ~50000 time units. The dashed line has a slope -1/3. 

hand much more pronounced than for the usual "GOY" model. This is related to the fact that 
the motion for the present model is also strongly intermittent and leads to corrections to the 
Kolmogorov theory. A study of the intermittency effects is reserved to a forthcoming publication. 

4 The theoretical predictions by Ruiz and Nelson. 

Ruiz and Nelson S, |j| have proposed a theory for the dependence of the mixing time on the 
hydro dynamical parameters. In the case of the passive scalar, a = 0, there is not a peak in the 
spectrum but an inhomogeneity created at the large scale will still persist for a mixing time which 
is composed of three terms 

Tpass=to + Hkd/k d ) + jj^y ■ ( 24 ) 

Here, the first term is the time it takes for a perturbation created at the large scale to reach the 
dissipative wavenumber kd (the Kolmogorov length). The second term is the time it takes for 
a perturbation to go from the dissipative wavenumber, through the "viscous-convective" regime 
down to the Batchelor wavenumber k' d = k^Pr 1 / 2 , and the last term is the time actually needed to 
dissipate the disturbance at k' d . The second term is found to be of the order 2 ^/ 2 ■ For large values 
of Re, the mixing time is therefore of the order to, unless InPr ^> Re 1 / 2 is which case the second 
term dominates. In the case of an active scalar on the other hand, a ^ 0, a peak occurs in the 
spectrum at a specific wavenumber k*, and the last term will therefore dominate when Pr 3> Re. 
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The corresponding mixing time is || 

1 Pr 

Active - ~ to— (25) 

Also, the wavenumber of the peak is predicted to be located at k* ~ koRe, where ko is the 
wavenumber of an initial perturbation. We observe from this theory, that for large values of the 
Prandtl number, the mixing time is much longer in the active case than in the passive case. 

We have tested the prediction of the theory given by Eq.([25|) using the shell model introduced 
in the previous sections. Fig. || shows a plot of the mixing time r ac ti ve versus the Prandtl number 




1 2 3 4 5 6 

Logio Re 



Figure 5: The mixing time T ac ti ve versus the Prandtl number Pr estimated from numerical calcu- 
lations of shell model with the following parameters: Upper curve: ./V = 16, ko = 2 -4 , v = 10~ 4 , a 
= 1CT 8 and the gate ip G = 1(T 6 . Lower curve: N = 19, k = 2" 4 , v =1CT 6 , a = 1(T 12 , and the 
gate ipa = 10~ 4 . The slope of the dashed line is 0.9. 

for two different values of the Reynolds number. In the first case, Re ~ 10 6 and Pr in the range 
from 10 2 to 10 6 and in the second, Re ~ 10 4 and Pr in the range 10 to 10 5 . The mixing time is 
estimated in the following way. At the shell corresponding to the k n -v&lue on which the maximum 
of the peak is localized, we monitor the value of <| ip n |>, where the average is over a specific time 
interval. As this value decreases below a chosen gate ipo, the corresponding time is associated with 
the mixing time T ac u ve , at that particular value of the Prandtl number. Then the value of Pr is 
changed and starting from the same initial conditions, using the same value of the gate, one obtains 
the new value of T ac u ve , and so on. The mixing time clearly diverges with the Prandtl number as 
a power law, T ac n ve ~ Pr 13 . The best fit to the data produces a value f3 ~ 0.9 for Re ~ 10 6 (as 
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indicated by the dashed line) and (3 ~ 0.95 for Re ~ 10 4 . The prediction (^) of Ruiz and Nelson is 
valid in the limit Pr 3> Re and this only holds for our data where Re ~ 10 4 , so the agreement with 
the theory is reasonable. It is however tempting to conjecture that the presence of intermittency 
might cause the exponent (5 to decrease below 1. The shell model presented here in Eqs. (13) 
and (|l^) exhibits strongly intermittent motion where the laminar periods are interrupted by bursts 
of violent motion. The higher the Reynolds number, the more pronounced is the intermittency 
and since the largest deviations from the prediction ( p5| ) is observed at Re ~ 10 6 , the effects of 
intermittency could diminish the value of mixing time due to the presence of the long laminar 
periods. 

Goldburg and co-workers investigated in several experiments the possible existence of the active 
coupling term in binary mixtures. In Ref. Q the presence of an active term was expected because 
the mixing time was found experimentally to increase dramatically with the value of the Reynolds 
number. Nevertheless, the experimental data did not follow the prediction (|25"|), Later experiments 
showed, however, that the long mixing times were caused by the fact that stirring of the binary 
mixture cools it down below the critical temperature, into the region where the system is immis- 
cible and phase separation is favored |J. Subsequent measurements on phase separations and the 
correlation functions of temporal fluctuations did not show any sign of an anomalous peak in the 
spectra @, §. In Ref. § it was already argued that intermittency effects might strongly influence 
the critical fluctuations leading to correlation functions which are stretched exponentials. That 
intermittency effects are very important close to the critical point is in accordance with the results 
presented in this paper. Conclusively, one must therefore say that, in spite of several experimental 
attempts, there is no clear evidence for the active coupling term in We return to this point in 
sections 7 and 8. 



5 The continuous cascade model for hydrodynamics 

The Kolmogorov scaling behavior is a static solution of the energy cascade model. In this section 
we discuss time-dependent generalizations of the simple A; -5 / 3 Kolmogorov behavior of the power 
spectrum, using a continuous version of the cascade model first discussed by Parisi |2Sfl . Initially 
these spectra start out as an arbitrary power behaviour A; 2p_1 , where p is some constant which can 
be selected as one wishes (if one includes necessary cutoffs in /c-space), but it turns out that after 
a short time the large k behavior becomes of the Kolmogorov type, with some time dependence. 
Thus, at a given time the spectra have a "two-slope" structure with a smooth interpolation. The 
solution is only valid in the inertial range, where diffusion can be ignored. 

The simplest hydrodynamical cascade model is given by the equation of motion |2C| , 

S 1 s 

(du n /dt + vk 2 u n )* = -ik n (u n+ iu n+ 2 u n _iu n+ i ^— u n -iu n - 2 ) + F n , (26) 

where F n represents an external forcing. From the point of view of energy (i.e. \ u n\ 2 ) conservation 
the parameter 5 is arbitrary, but it can be fixed by requiring conservation of (generalized) helicity 
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|23| |. In our case we keep the parameter 5 arbitrary. 

Some years ago Parisi [28] studied eq. (26) in the limit [] r — ► 1, meaning that the distance 



between the shells goes to zero. Taking r = 1 + e it is easily seen that one gets 

I- + uk 2 ] u* = -ie(2 - 5)k (u 2 + Sku^f] + F(k) + 0(e 2 ). (27) 
at J \ ok J 

One might consider this equation to be a "model of a model", and thus very academic. However, as 
we shall see, this model and, in particular, its generalizations, satisfy the same conservation laws as 
the corresponding discrete models. Therefore one can equally well consider the continuous version 
as a model in its own right, with the advantage that it is considerably simpler than the discrete 
versions. 

We can now scale e(2 — 5) into time t and then let e — > 0. As just mentioned, the resulting 
model can be considered as being independent of the discrete version, since it satisfies the relevant 
conservation laws. Instead, one can consider eq. fl27| ) to be an approximation to the discrete model, 
to be supplemented by higher order terms in e if needed. 



We shall study eq.(|27[) in the inertial range, where viscosity can be ignoredfj. Also, we disregard 
the possible forcing term, so it is then clear that the motion must die out after some time if diffusion 
is included. Hence we study the equation 

^ = -ik(u 2 + 3ku^)+0(e). (28) 



dt V dk J 

We now choose a special phase and make the ansatz u = ik p f(k q t). Inserting this in eq. (|28| ) we 
obtain q = 1 + p, i.e. 

u = ik p f(k 1+p t). (29) 



This scaling was first considered by Parisi |28|| , except for the special case p=l, which was introduced 
many years ago by Heisenberg [pCj] . It solves the discrete as well as the continuous model. 

A scaling of the type (^) is consistent with the well known invariance of the Navier-Stokes 



equations (see for instance |31[|), 

I — > XI, \ h u, t -» X 1 ~ h t, and v -> \ l+h v, (30) 

in the inertial range, where we can take u=0. The reason for this is that the self-similarity ( |30| ) 
can be translated to fc-space with I — > 1/k, and it then corresponds to the scaling ([29] ) with 
p = —h. The main point is that the scaling variable k 1+p t is then an invariant. Also, in the scaling 
(|30|) I is usually interpreted as the scale of an eddy, and the typical velocity of this eddy is then 
|u(x + 1) — u(x)|. This compares excellently to the velocity mode u n used in A:-space, where u n 

1 In ref. Q this limit was considered for a different type of shell model. The resulting equation is linear, in contrast 
to those considered in this paper. 

2 There exists the possibility of taking the more exotic limit |e(2 — S)\ — > a finite value, so that |<5| approaches 
infinity. 

3 For the special case p—1 we shall include diffusion later in this section. 
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is the velocity increment over an eddy of scale I ~ l/k n . Let us further note that for p = —h=l, 



diffusion can be included, since (30) then leaves v invariant. 



It should be emphasized that the expression (^9|) has the explicit power k p in front to accom- 



modate the self-similarity transformation (30) for u. However, the function / depends only on the 
quantity k 1+p t, which is invariant under the self-similarity transformations (|30|). In the absence of 
a solution for /, this function can be completely arbitrary from the point of view of self-similarity, 
and hence e.g. the time evolution cannot be predicted at all. This implies that a priori there is no 
agreement with K41 theory (see [31] for a general discussion of K41). Of course, having a solution 
for / changes the situation. 

Eq. (|29| ) thus means that the velocity mode is initially assumed to be of the form k p . Physically, 
one can imagine that this initial condition is produced by some external force. In this sense the 
selection of initial conditions is equivalent to initial forcing. Also, p governs the initial correlation 
function < Ui(x.)uk(y) >. The cases p=3/2 or 1 correspond to Gaussian disorder in three and two 
dimensions, respectively, i.e. < itj(x)«fc(y) >oc <5jfc(5 3 (x — y) or <5jfc<5 2 (x — y), respectively^} This is 
because the k-space energy spectrum is given by 

E(k,t) = \u(k,t)\ 2 /k = const k D ~ l J d D xexp(ikx) < u(x, t)u(0, i) >, (31) 

in D dimensions. 

In order to have a convergent energy 



E = J dk\u\ 2 /k, (32) 

we obviously need a ultraviolet cutoff for p >l/2. Similarly, for p <0 an infrared cutoff is needed. 
At this stage we need a discussion of the boundary conditions associated with (|27|) and (|28|). 



From (27) we get energy conservation in the absence of forcing and viscosity provided 

k\u(k,t)\ 3 -> for k^O and oo. (33) 

If this condition^] is not satisfied for k — > oo there is "diffusion at infinity" . 
Inserting (|29|) in eq.(p3), we obtain the following equation for /, 



df(x) (1 + 3p)/(x) 



dx 1 + 3(1 + p)xf{x) ' 

This equation can be simplified by the substitution 



k 1+p t. (34) 



f(x) = g(x)/x, (35) 



4 Here we leave out a discussion of the consequences of divu=0, which requires a projection operator in the 
definition of Gaussian randomness. This is of no relevance in the following. 

5 In the discrete version ( ]26| ) there is, strictly speaking, a similar boundary condition if n goes to infinity. This 
is because when one checks energy conservation, sums of the type k n uuu (the u's have different indices) are 
encountered. Although there is a complete cancellation of these terms, the sums only exist in a strict mathematical 
sense if all terms of the type k n uuu vanish for n — » oo. 
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and we get 

dg{x) g(x) + 2g{xf 



(36) 



dlnx l + 3(l+p)g(x) 
The solution is given by 

g(x)(l+2g(x))^y 2 =x/x . (37) 
Here xq is an arbitrary constant which gives the strength of the initial velocity mode, 

u(k,0) = ik p /x . (38) 

When x — > oo we find from ([37]) that 

g ( x ) _> 2 -(l+3 P )/3(l+ P ) ( x/xo )2/3(l+ P )_ (39) 

Inserting this in ((2^) and ([??]) we get 

_ l+3p_ , (l+3p) i 

|li(M)l -»2 Xo 3 ^+W t 3(1+,,) (40) 

Thus we see that irrespective of the initial spectrum ( p8[ ) the velocity approaches the Kolmogorov 
spectrum with a time dependent amplitude for large values of k and/or time. Note that this decay 
law does not agree with the classical theory put forward by Karman, Howarth and Kolmogorov 
which is of course due to the non-triviality of the function / in (^) ||32"| , [l0|| (see the general 
discussion on classical results on decay laws by Frisch |MJ ) . 



There are a few special cases where eq.(37) can be solved explicitly. The simplest is the case 
where p = —1/3, where we get g(x) = x/xq, leading to the time independent Kolmogorov spectrum, 



u 



l t)\ = k-V 3 /x , for p = -l/3. (41) 



This result is trivial, since it is easy to see that the original equation ( |28| ) has ( |4l| ) as a static 
solution. 

A non-trivial result can be obtained by considering the case p = 1/3, where eq.(|37|) becomes 



second order in g. Using (29) and ( |35| ) we then obtain 

KM) I = ^ ("I + yJl + W/H/xo^) . (42) 

For small but non- vanishing t, the slope in the corresponding power spectrum changes from -1/3 
to -5/3. For k — > oo the boundary condition (^) is not satisfied. This is simply a consequence of 
the fact that the Kolmogorov spectrum does not satisfy this condition. As mentioned before, we 
only expect the solution ( f42|) to be relevant for the inertial range, where diffusion is negligible. For 
completeness we mention that the energy integrated to h = oo is given by 

E = V2/(x 2 Vt). (43) 



The energy contents approaches zero as time goes on simply because the class of solutions ( 29 ) 
correspond to an inverse cascade moving towards smaller values of k, ultimately reaching the value 
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except in the point k = 0, where it is given by 1/x 2 ,. In "x-space" this means that small scale 
structures in time become structures of infinite extension, with no local energy. This shows up if 
we consider the "integral scale" Iq used in turbulence theory as a measure of possible large scale 
structures, 

l = J ^-E{k) (J dkE(k)*j , with E(k) = \u\ 2 /k. (44) 

This quantity behaves like t 3 / 4 , and in the general case it goes as fVCl+P). Hence, for large times 
the structures become very extended. 

Another case which can be solved is p = 1, since eq.(|37|) then becomes a cubic equation for g. 
We find 

\u(k,t)\ = R(s)/kt, (45) 



R ( s ) = l\ [ 1 + s [ 1 + \l l + 1 ) \ + + + r ) i " 2 I> ( 46 ) 



where 

RCsl = 

6 
and 

„2 




s = 27k z t/x . (47) 

The energy spectrum initially has the slope +1. For small k the slope remains +1, but at larger 
values of k the slope turns into -5/3, in accordance with the universality of the Kolmogorov spectrum 
at large fc's. Again there is an inverse cascade moving the energy towards smaller fc's as time passes. 

The p=l case is interesting from the point of view of studying the effect of diffusion. The main 
feature is that in general the ansatz (|2^) is not consistent with the viscosity term on the left hand 
side of eq. (p7|). However, when p=l, the powers of k nicely divide out on both sides of eq.(p7|), 
leaving them as functions of the scaling variable x only. Instead of eq. ( p4| ) we obtain 

df(x) = 4/(x) 2 + uf{x) 
dx l + 6xf(x) ' 1 ' 

The substitution (|35|) cannot be used to solve this equation because of the dissipative term. How- 
ever, eq.([48|) can of course easily be solved numerically, and the result compared to the analytic 
solution (|46|). To see the effects expected, let us compute the diffusion cutoff defined by f 2 ~ vf 



(see eq.(|48D), using eqs. (gSj) , (|35[) , and (|39| ) applied to the case p=l, 

x D ~ x Q 1/2 z," 3 / 2 , or k D ~ ^ 3 / 4 x- 1/4 t" 1 / 2 , (49) 

which is what one one expects as far as the dependence on v is concerned. However, the time 
dependence should be noticed. The latter reflects the fact that we have an inverse cascade. The 
cutoff (^) is in qualitative agreement with what one obtains by comparing the analytic solution 
(fl6|) with the numerical solution of (|48|) . In fig. ^|we give an example of a solution, where we plot 
the energy. It is clearly seen that there is a change from the initial slope +1 to -5/3. 



For a general p the result (49) can be extended to 

k D ~ v^l A xo 1/2(1+P) t-d+3p)/4(i+ P ). (50) 
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Figure 6: A double logarithmic (base 10) plot of the energy spectrum E(k, t) as a function of k in 
the continuous model. In the left hand figure we show ten small values of time, from t=0 (where 
the E goes like k) to t=10~ 2 . In the right hand figure we show E for ten values of time between 
t=l and t=l.l. In both cases time moves from the right to the left. The viscosity is v=l0~ 2 , xq 
=1, and from (^) logio^D is of order 2.5 (left) and 2 (right), in reasonable agreement with the 
figures. 



The behavior of ko as a function of v is thus universal, i.e. independent of p, whereas the time 
dependence of kp is clearly non- universal. 

The results (|49| ) and (|50| ) are to some extent consistent with Kolmogorov's scaling arguments 
as far as the ^-dependence and as far as the time dependence are concerned: The energy should 
go as E(k) ~ e 2 / 3 /c~ 5 / 3 , where e is the kinetic energy per unit mass (anyhow put equal to one) 
and time, with e = dE/dt ~ ku 3 for large k. Now from (f40|) ku 3 ~ ^( 1 + 3 p)/( 1 +p) f or large 
k. From Kolmogorov's arguments one expects the dissipation scale ku k[) ~ 

^-3/4 ^-(l+3p)/4(l+p)_ ^^jg i s exactly the result (pC|), also as far as the time dependence is concerned. 
Of course, the fact that e is not a constant in our case, does not conform to Kolmogorov's argument. 

In the case p =1 it is possible to obtain an implicit equation for g(x) even when diffusion is 
included. For v 7^0, proceeding like in eqs . (|34|) , (|35|) , (|36|) , (|37|) , and (|46|), we get that eq. (|37|) is 
replaced by, 

kH ( r k2t 1 \ 
g(k 2 t)(l + 2g(k 2 t)f = — exp \-v / dx . (51) 

xo \ Jo l + 2g{x) J 

Therefore eq. (|4^) is changed to 

\u(k,t)\ = R(d)/kt, (52) 

where R(d) is given by (f46"[), and where 

x \ Jo 1 + 2^(x) y 
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This equation can be used iteratively, starting from the behaviour for i^=0. To the lowest non-trivial 
order in v eq. ( |52|) becomes 

|u(M)l = R (27k 2 texp(-vk 2 t)/x ^ /kt. (54) 

for k 2 t small. It is interesting that ( |54} ) is a self-consistent solution of eq. ( ]5l| ) for fc 2 i small and 
large. Thus the decay of the velocity at large k is exponential. For k 2 t neither small or large, the 
expression given by (|46| ) and (|53| ) indicates that an exact inclusion of diffusion is rather complex, 
even in this simple model. 

It is very interesting that the scaling behaviour for the case p= 1, i.e. u(k,t) = kf(k 2 t) was first 
considered by Heisenberg in his model, where in the range from to k, the action of all smaller 
eddies are assumed to be represented by an effective viscosity [j30| l. Although it is not obvious 
that the cascade model satisfies this assumption, the Heisenberg scaling appears as a solution. 
Physically the fc 2 £-scaling can be understood |3(J by assuming that the spectrum is determined by 
one length only, namely the length ~ 1/ko of the largest eddies. Let the velocity of these eddies be 



vq. On dimensional grounds it then follows that [30] 



— (1/ko) ~ v , and ^(V u o) ~ h, (55) 



from which one gets ko ~ v$ ~ 1/yt. This is exactly the scaling in the p=l case. However, the 
continuous model has an exponential decay for large k, in contrast to Heisenberg's model, which 
has a power behaviour, presumed to be unrealistic |[^]. 

The solutions discussed above refer to the case when the phase is fixed. One can try to find 
solutions where the phase plays a dynamical role by making the scaling ansatz u = ik p f (k 1JrP ) with 
/ a complex function, in which case eq. (^4[) is replaced by 

^M. = _ 3{1+p)xfix) Vp_ _ (i + 3 p)f(x) 2 . (56) 
ax ax 

For the case p=l/3, introducing again the substitution (35) with g(x) complex, this leads to the 

result 

g*(x) + 2g(x) 2 = x/xq, x = k 4/3 t. (57) 

This equation has the previously discussed real solution corresponding to (|42|), as well as a new 
complex solution with 

Reg(x) = l/A, lmg(x) = J (3 - 8x/x Q )/4. (58) 

This solution is obviously only valid for x < 3/8 x$. Since x can take any value from to oo the 
complex solution should be rejected. 

To summarize the results obtained so far, one can say that the continuous model gives interest- 
ing, non-trivial results. This model is most interesting for the case where p=l, since then diffusion 
can be included. The initial condition \u(k, 0)| oc k then corresponds to an initial Gaussian dis- 
order in two dimensions. Therefore we believe that this model may be of most relevance in two 
dimensions. This is consistent with the fact that the model has an inverse cascade. 
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There exists a generalization to a model with helicity, namely the continuous version of a model 



introduced by Biferale and Kerr [29|, leading to |C 



d_ 

dt 



+ vk A {u + y 



-ik I 4ku~^r- + 2ku + ^- + (2 + a)u + u~ - a(u~) 2 



dk 



dk 



(59) 



There is a similar equation with +<->-. The energy and the generalized helicity are conserved, 



E = J ^(|u+| 2 + |tT| 2 ), H = J ^k a {\u + \ 2 - |n"| 2 ). 
Making the scaling ansatz 

u+ = ik p f(k 1+p t) and vT = ik p h(k 1+p t), 

one obtains from (^) by ignoring viscosity 

df(x) _ 2(1 + p)xfdh/dx + (6p + 2 + a)fh - ah 2 
dx \ + A(\+p)xh 

There is a similar equation with / and g interchanged. These equations have a potentially much 
richer structure than eq.(|34|). We also mention that the discrete and continuous GOY equations 
have been generalized to magnetohydrodynamics, using a helicity decomposition [34]. 



(60) 



(61) 



(62) 



6 The continuous shell model for turbulent mixtures 

In section 2 we discussed the shell model for binary mixtures. Proceeding exactly as in the previous 
section, we can now derive the continuous version of the relevant equations. Here we shall just give 
the results. For convenience we define <f> = kip, and after a rescaling of time the equations become 

% + D*V = % k {-v4> + #g + 2vk^ , (63) 

and 

+ v k 2 v* = ik (2a(/) 2 + c(v 2 + Zvk^J ^ . (64) 
Here v is the velocity mode and c is an arbitrary constant. These equations conserve 

E - / f h/f , and £(k| 2 + a\M 2 ) - J f (M 2 + ^\ 2 ), (65) 

provided we have the boundary conditions 

v\<p\ 2 /k -> 0, k\v\ 3 -> 0, k\<j)\ 2 v -> 0, (66) 

for k — > oo and 0. Again there may be "diffusion at infinity", as discussed in the last section. 

We can now look for scaling solutions of eqs.((6^) and (|64]) in the inertial range. It turns out 
that they should have the form (analogously to MHD) 

(f>(k,t) = ik a 4>{k l+a t), and v(k,t) = -ik a v(k 1+a t). (67) 
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These results have been found by requiring that powers of k on the two sides of eqs. (63) and (| 
should cancel out. This fixes the powers in v and (j) to be identical. It is interesting that these results 
are again consistent with the exact self-similarity of the Navier-Stokes eq. (|30|) , supplemented by 
ip — > X~ a ip, as can be seen from eq. (pi). From eqs.(|63[) and (p3) we then get 



_ (f)'(xY = (3a - l)v(x)(j)(x) + (1 + a)x(f>(x)v {x) + 2(1 + a)xv(x)<j/(x) (68) 

and 

v'(x)* = -2a4>(x) 2 - c((l + 3a)v(x) 2 + 3(1 + a)xv(x)v'(x)). (69) 

Here the scaling variable x is given by k 1+p t. 

In general the scaling ansatz ([67]) is inconsistent with diffusion, and hence eqs. ([38]) and ( |69| ) 
can only be used in the inertial range. However, as already seen in the last section, the case a=l 
is an exception. In this case all powers of k neatly cancel out even in the presence of diffusion. If 
we assume that the functions <f) and v are real, the equations become 

- 4>\x)/4>{x) = (2v(x) + 2xv{x) + D)/(l + 4xv(x)), (70) 

and 

v'(x) = -(2a<j)(x) 2 + Acv(x) 2 + vv(x))/(l + 6cxv(x)), (71) 

where the scaling variable is now x = k 2 t. 

Eqs. ( |70D and ( |7l| ) can be reformulated in a way which is similar to what was done in the last 
section. From ( [70[) we get 

* [kH) = 7TTwmm>r v (r°f dx/(1 +4I5(I)) ) ' <72) 

Similarly eq.(^) can be integrated to give the implicit equation 

+ ^t^t)? - C (0)exp (-^fc—L— + , (73 ) 

provided v(0) does not vanish. This result is similar to eq.(^), and a comparison shows that the 
last factor on the right above (containing cp 2 /v) is very similar to diffusion. Thus in this particular 
model, the effect of (f> on the velocity field is essentially to provide some additional diffusion. 
Comparing with the last section, we therefore expect that the Kolmogorov regime for the velocity 
mode kv becomes smaller, since the effective diffusion becomes stronger. 

In the inertial range eqs. ( |72|) and (|73| ) can be explicitly solved analogously to what we did in 
the last section. Eq. ([73]) has the solution 

v{k 2 t) = R(27cv{0)k 2 t)/ck 2 t, (74) 

where R is given by eq. (|46|). Eq. (|f^) then gives 



(j)(k 2 t) = 0(0)/ Jl + 4R(27cv( y 0)k 2 t)/c. (75) 
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Like in the last section it is easy to see that for large k 2 t eqs. fl74| ) and (75) imply the Kolmogorov 
and Obukhov-Corrsin exponents for the velocity and for the field <j) = ip, respectively. 

The advantage of the basic equations (^) and (|7l|) is that they are just two ordinary coupled 
differential equations. Therefore, they are much simpler than the usual large number of coupled 
cascade equations. However, it should be remembered that the set-up is very special. We need to 
assume initial (t=0) spectra which are linear in k. Thus, it is not possible to start e.g. with an 
initial Gaussian spectrum. However, the results in the discrete cascade model reported in section 
3 indicate that the final results are independent of the initial state. However, in any case one 
may wonder whether the simplicity has not been achieved at the cost of loosing the physics of the 
problem. This will be discussed in the next section. 

7 Results in the continuous model 

In this section we shall compare the continuous cascade model with the discrete one by obtaining 
qualitative and quantitative results. From the scaling which we introduced for v and (j), it is clear 
that the continuous model will produce an inverse cascade, whereby energy is transferred from 
large to small fc-values. Such a phenomenon is known in magnetohydrodynamics (MHD) for the 



magnetic energy [26] and in an MHD discrete cascade model p7[| , where it is presumably due to an 
inverse cascade in the three-dimensional magnetic helicity (for another shell- model on this point, 
see J33|). The continuous cascade model for MHD gives results very similar to the discrete 
model. In this connection, it is of interest that the basic equations used for turbulent mixtures are 
rather analogous to the MHD-equations. 

We begin by a qualitative discussion of the results expected, based on eqs. and (|73|). If 
diffusion is ignored (and a = 0) these equations are explicitly solvable. Eq. ( |74| ) implies that v is 
given exactly like the velocity u in section 5, so 

xv(x) oc x 1//3 (76) 

for large values of x = k 2 t in a range, where diffusion can still be ignored. From eq. (f75|) we then 
see that the quantity C(k) = ip 2 /k = (j? /k behaves like /c _5//3 with the Obukhov-Corrsin exponent 
(with some time dependence). If the Prandtl number Pr = v/D is large, diffusion in the first 
place acts on the velocity. Thus, xv starts to decrease, and from ( |75|) it follows that (j) approaches 
its initial value, i. e. C(k) goes like 1/k (with some time dependence), which is the Batchelor 



behaviour. Therefore, without doing any numerical calculations, we see from eqs. (73) and (75) 
that C(k) must change slope from -5/3 to -1. When k becomes so large that the diffusion governed 
by D is operative, of course C(k) decays exponentially, according to eq. ([7^). If, on the other hand, 
Pr is of order one or less, it follows by the same reasoning as given above (based on eqs. ([72]) and 
(f73|)) that only the -5/3 slope materializes itself, before the exponential decays set in. This is in 
agreement with the results obtained in the discrete model, as discussed in section 3. 
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We have made some numerical calculations, using the values 

v = 1CT 1 , D = 1(T 7 , a = u 2 , and Pr = 10 6 . (77) 

The initial values are given by u(0)=l and 0(0)=. 0001. We also took the constant c in eq. ([n]) to 
be 1/2. The results are presented for the spectrum C{k) = (f> 2 jk = if) 2 jk and the velocity mode 
kv. In Fig. [7] we show C(k) for relatively low and for relatively large times. The initial behaviour 




Figure 7: A double logarithmic (base 10) plot of the spectrum C{k) as a function of k in the 
continuous model for (a) ten small values of time, from t=0 (where the C goes like 1/k) to t=10 -2 , 
and (b) for hundred times larger. Time moves from the right to the left. 

is linear in k (from our initial condition that 0(0) is a constant). From Fig. ^ we see this behaviour 
for i=0. When t increases, the slope changes to -5/3 in a range of k- values from approximately 
10 to slightly less than 10 2 . From Fig. || for the velocity mode, we see that at a k- value around 
100, the velocity decreases rapidly. Therefore, as one can see from Eq. (f72|) and as one can also 
see from Fig. 0(a), after the initial decrease with slope -5/3, the spectrum C(k) increases until it 
regains the Batchelor slope -1. Ultimately, for k slightly above 10 4 one sees that C{k) starts to 
decrease exponentially. This repeats itself at later times, as is also seen from Fig. |7[ but when time 
increases, the "velocity of change" k/t oc \j\ft decreases, so the curves for the ten different times 
in Fig. ^ are much closer at later times. In the region where C increases (between slope -5/3 and 
slope -1), there are very short range correlations in x-space. This could be a rudimentary version 
of an intermittency fluctuation. 

We have also investigated the question of equipartition. For small k the kinetic energy domi- 
nates. Around the value of k where the peak appears, there is equipartition in essentially only one 
point. The energy E^ = akif) 2 then dominates for larger values, where the kinetic energy becomes 
very small compared to E^. It should also be noticed that the peak is time dependent, and moves 
towards smaller values of k with a "velocity" k/t = \j\ft- 
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Figure 8: Double logarithmic plot (base 10) for the velocity mode v(k,t) = kv as function of k for 
(a) ten time values between (where v{k, 0) = k) and 10 -2 , and (b) for 100 times larger times. 



Finally, we have studied the influence of the parameter a. It seems that the continuous model 
differ from the discrete one with respect to this point. For the various values of Pr and v we have 
studied, we find that the influence of a is only rather marginal, and we do not find any really 
spectacular effect of a. It would therefore be interesting to study the case where the initial velocity 
is zero (or very small). In this case there is no effect (or only a slight effect) if we have a passive 
scalar, so with a non-vanishing one will see an effect (almost) entirely due to a. This problem can 
be analyzed from eqs. (|72] ) and (ff^) by replacing v(0) by v(e) and performing the limit e — >0 in 
such a way that v approaches zero. We can also study this by analyzing eq. ( |71~D near x=0, using 
the boundary condition -u(0)=0. We get from |n]) (compare with eq. (|l2|)) 

v(x) « -2q0(O) 2 x + 0(x 2 ), (78) 

Thus the velocity scaling function must be negative and linear for small values of the scaling variable 
x = k 2 t. 

We have also studied this problem numerically. In Fig. ^ we show the resulting scaling function 
v 2 = v(k,t) 2 /k 2 as a function of the scaling variable. The values of v and Pr are as used before, 
but initially we assume u(0)=0, and (f)(0) = ip(0)=.0005. Although the velocity-scaling function is 
very small, there is clearly an effect. However, the resulting back-reaction on <j) is small in most 
cases. This can be seen from eq. ([72]), since the velocity field in the denominator on the right hand 
side is small. However, there may exist initial conditions where this is not true, and where xv(x) 
may approach -1, causing tjj to diverge. In such a case stabilizing terms of higher orders in t/j must 
be included in the basic eqs. (||) and @. Disregarding this possibility, the density function ip is 
rather insensitive to velocity fluctuations, when the kinematic viscosity v is much larger than the 
viscosity D of ip. This effect is somewhat similar to what has recently been seen in MHD |34|, |35| for 
large "Prandtl" numbers (i.e. for the kinematic viscosity much larger than the Ohmic diffusion). 
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Figure 9: A log-log plot (base 10) of the scaling function v 2 = (v(k,t)/k) 2 as a function of the 
scaling variable x = k 2 t for the initial conditions v(0)=0 and <^>(0) = ?/>(0)=.0005. 



From an experimental point of view it would be interesting to study the case where the initial 
velocity vanishes or is small, and where some initial distribution (e.g. random) is established with a 
non-constant gradient of tp. The fluids should then be set "spontaneously" in motion. The gradient 
field is analogous to the magnetic field in MHD. Therefore, such an experimental set up would be 



somewhat analogous to the study of primordial magnetic fields in the early universe (see refs. [27] 
and p^l), where the magnetic field can induce a velocity field. It would clearly be of interest to 
investigate this analogy in an earthbound laboratory. 

To conclude this section, we do not find the large peak predicted in ref. || for increasing Prandtl 
numbers. There is a small time dependent peak of magnitude less than a half decade, as can be 
seen e.g. from Fig. ^. This is to be contrasted with a peak of several decades in ref. ||. Whether 
this is a shortcoming of the continuous model is, of course, an experimental question. However, our 
results are not so different from those of the discrete GOY model, discussed in section 3, where the 
peak is of the same order of magnitude as found here. However, the nature of the peak is different 
in the two cases, since the a dependence differ, as mentioned above. 

8 Conclusions 

The main result of the present paper is that intermittency effects are likely to play an important 
role in turbulent binary fluids. The influence of intermittency is not small; compared to previous 
studies of non-intermittent binary fluids || [| the peak in the concentration spectrum is much 
less pronounced and less persistent because the fluctuations tend to "surround" and diminish the 
peak. It is therefore not easy to see the effect of the active coupling term (||) from the spectra, 
in accordance with experiments j§-[H]. The transport coefficient a in (13) presumably has its most 
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dramatic effect in the case where the initial velocity vanishes. Here the existence of a non- vanishing 
a implies the "spontaneous" generation of a velocity field, provided there is an initial variation in 
the gradient of ifi. The observation of such an effect would have an analogy in MHD, where an 
initial ("primordial") magnetic field induces a velocity field, which may be of relevance in the early 
universe. On the other hand, if such an effect is not observed in binary mixtures, this would 
indicate that the active coupling term (||) is probably not present. If so, binary mixtures would 
not be analogous to MHD. 

Also, we see that the continuous model gives results which are rather similar to the discrete 
model. In the continuous case there also exists a time-dependent peak. 

We thank David R. Nelson and Walter Goldburg for stimulating discussions. We also thank 



George Savvidy for reminding us of Heisenberg's paper [30], and Yuri Makeenko for much help in 
inserting the figures. 
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